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Abstract. As a result of the application of coarse-graining procedures to describe 
complex fluids, the study of systems consisting of particles interacting through 
bounded, repulsive pair potentials has become of increasing interest in the last years. 
A well known example is the so-called Generalized Exponential Model (GEM-m), 
for which the interaction between particles is described by the potential v(r) = 
eexp[— (r/c) 1 ™]. Interactions with m > 2 lead to the formation of a novel phase 
of soft matter consisting of cluster crystals. Recent studies on the phase behavior 
of binary mixtures of GEM-m particles have provided evidence for the formation of 
novel kinds of alloys, depending on the cross interactions between the two species. 
This work aims to study the dynamic behavior of such binary mixtures by means of 
extensive molecular dynamics simulations, and in particular to investigate the effect 
of the addition of non-clustering particles on the dynamic scenario of one-component 
cluster crystals. Analogies and differences with the one-component case are revealed 
and discussed by analyzing self- and collective dynamic correlators. 
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1. Introduction 

Coarse-graining procedures provide a general way to represent mesoscopic aggregates of 
several architectures (e.g., linear and hyperbranched polymers, micelles, microgels, etc.) 
as particles interacting via an effective, isotropic potential V e s(r). When the center of 
mass of such aggregates can coincide without violating excluded volume interactions, 
V c ff(r) turns out to be ultrasoft and bounded [TJ |2j [3]. These characteristics may 
lead to cluster formation providing V e s(r) belongs to the so-called Q ± -class, i.e. its 
Fourier transform (FT) V c g(q) displays an oscillatory decay around zero (potentials 
with positive FT are referred as Q + potentials) [H [3] . Recent numerical investigations 
of amphiphilic dendrimers and of ring polymers with different degree of knotedness 
demonstrate that such macromolecules are candidates for experimental realizations of 
the former clustering scenario, since their effective interactions are indeed of the Q ± -class 
El [8] . A particular theoretical realization of Q ± potentials is given by the generalized 
exponential model (GEM-m), for which the interaction is described by the potential 



where r is the interparticle distance, a is a measure of the particle size, and e is the 
interaction strength. As density increases, the GEM-m system shows, for m > 2, a first- 
order transition from a fluid of clusters to a cluster bcc-crystal above the triple point 
temperature, and a subsequent structural phase transition from a cluster bcc to a cluster 
fee crystal. This occurs at arbitrarily high temperatures [91 [10], contrary to the case 
of Q + -systems, which exhibit reentrant crystallization in the temperature-density plane 
[TT] . A remarkable feature of the cluster fee crystals is that upon increasing the density 
the lattice constant approaches a constant value af cc = 2n\^3/q*, with q* the wave vector 
at which v(q) takes its absolute minimum. An obvious consequence of the density- 
independent value of af cc is that the population of each lattice site scales proportionally 
with density. These properties are rather different from usual crystallization features 
in colloidal systems. Moreover, these differences are not restricted to the structural 
properties but also to the dynamical ones. Regarding the slow dynamics in such cluster 
crystals, incessant hopping between clusters has been revealed, which fully changes the 
initial identity of the clusters without altering the lattice structure. Another peculiar 
dynamic feature of these systems is a decoupling between self- and collective out-of- 
lattice correlations [121 EE] • 

Recently, the phase behavior of mixtures containing a non-clustering component 
(m = 2) and a clustering component (m = 4) has been investigated within the framework 
of the Density Functional Theory (DFT) [T3J [15] . For the case of non-demixing systems, 
evidence was found for the formation of novel kinds of alloys, i.e. mixed cluster crystals 
[13] . Though static features of these mixtures have been studied in detail, to the best 
of our knowledge no information has been reported on the corresponding dynamical 
aspects. This article aims to shed light on such aspects. By means of extensive molecular 
dynamics (MD) simulations, the findings displayed by one-component cluster crystals 
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[T2~| ITS] are extended to mixtures of clustering and non-clustering particles. Analogies 
and differences with the one- component case are revealed and discussed by analyzing 
self- and collective dynamic correlators. The article is organized as follows: Model and 
simulation details are given in Section 2. Static and dynamic features are presented and 
discussed in Section 3. Conclusions are given in Section 4. 

2. Model and Simulation 

As a previous step, the stability limits of the homogeneous fluid phase (coexistence 
and A lines) were determined for representative examples of the former class of binary 
mixtures. These limits were obtained by solving the Ornstein-Zernike (OZ) equation 
within the mean field approximation (MFA) [16]. Based on such results, the specific 
systems to be simulated were determined and thus the parameters defining the mixture 
(relative particle size, mixture compositions, temperature and density, self- and cross- 
interactions) were chosen in such a way as to avoid the region of demixing in the phase 
diagram. Because of the high total densities to be considered, the former parameters 
were also selected with the goal of reducing computational expense. 

The investigated system was a mixture of big GEM-8 particles and small GEM-2 
(gaussian) particles. This mixture contains a non-clustering component (GEM-2, in 
the following referred as A-particles) and a cluster-forming component (GEM-8, in the 
following referred as S-particles). The dynamic features of this system were investigated 
at fixed temperature T for a broad range of densities p = pa + Pb an d several relative 
compositions x = Pb/p- The former densities are defined as p a = N a /V, N a being the 
number of particles of the species a and V the total volume. In what follows, we choose 
units in which Boltzmann's constant k-g has the value fcs = 1- The cross-interaction 
was chosen to be a GEM-4 potential, which also displays cluster formation. In this way, 
the following interaction parameters were adopted for the potential (JTJ: e a p = e = 1.0, 
o'aa = 0.3a, gab = 0.6a, <jbb = o = 1.0, and tuaa = 2, tuab = 4, tubb = 8. With 
these parameters, a common cutoff range R c = 1.5a was introduced in the simulation 
for all the interactions. 

MD simulations were performed in a cubic box with periodic boundary conditions. 
The equations of motion were integrated in the velocity Verlet scheme [17] , with a time 
step At/r ranging from 0.001 to 0.005, where r = y / mcr 2 /e and a common mass m — 1 
was used for all the particles. The size of the simulation box, L^ ox , was typically of 7 or 8 
times the lattice constant af cc of the considered mixture. The value of af cc was obtained 
from DFT as explained in [14], [15]. In the following, density, time, distance, and wave 
vector will be given, respectively, in units of a -3 , r, a, and a -1 . Initial configurations 
were generated by placing the particles uniformly at the sites of the fee lattice defined 
by af cc . An equilibration run was performed in which the system was thermalized at 
temperature T = 0.30 by periodic velocity rescaling. Typical equilibration times ranged 
from 2 x 10 5 to 2 x 10 6 time steps according to the studied density. After reaching 
equilibrium, manifested by the absence of any drift in internal energy and pressure, a 
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Figure 1. Left: A segment, at the high-x/moderate p region, of the phase diagram of 
the studied ultrasoft mixture. Empty and filled symbols indicate the simulated state 
points in the (p, x)-plane, at both sides of the melting line. The latter is displayed as 
a dashed guideline for the eyes. It is estimated from the emergence of Bragg peaks in 
the partial static structure factors (see below). The solid line is the A-line, obtained 
from the OZ-MFA equations. Right: Slabs of the simulation box for x = 0.65 and 
p = 3.50, parallel to the (OOl)-plane (top) and to the (lll)-plane (bottom). For the 
sake of clarity, only GEM-8 particles are displayed, in both frontal and lateral views, 
with a size much smaller than the real a. Particles in different layers are represented 
with different colors. 



production run was performed in the microcanonical ensemble, at different densities, for 
compositions x = 0.65, 0.80 and 0.95. Typical production runs ranged from 10 6 to 10 8 
time steps. 

3. Results and Discussion 

3.1. Phase diagram and static structure 

The partial structure factors S a p(q) (a, (3 G {A, B}) are related to the propensity of a 
fluid to sustain spontaneous density fluctuations of wave vector q. They are defined as 

( } = (p a (q,0)p/,(-q,0)) (2) 
jN a N p 

with p a (q, t) = X^iexp[iq • r"(t)] and the sum is performed over the coordinates r° 
of all particles belonging to the species ajij] According to the OZ equations they take, 

| For uniform phases, the statistical average (• • •) in the right hand side of eq. (|2|) above renders its left 
hand side a function of q = |q| only. For the crystalline phases to be considered in what follows, we have 
performed an additional rotational average over symmetry-related values of the reciprocal vectors, so 
that the resulting quantities (static and time-dependent correlation functions) are shown as functions 
of q only. 
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Table 1. For each investigated composition x, density p* for the A-line, parameters 
defining the fee lattice (D(q*) = 0, af cc = 2iry/3/q*, see [2]) and box size I/box used 
in the MD. 
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1.8960 
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0.95 
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5.84 
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1.00 
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5.86 


1.8571 


14.8571 



in the homogeneous phase, the generic form S a p(q) = N a p(q)/D(q) where, within the 
MFA, the denominator depends on the Fourier transform of the interaction potentials 
according to 

D(q) = 1 + ^U s (q) + {£)x{\ - x)U d (q), (3) 

where U s (q) = (1 - x)v AA (q) + xv BB {q) and U d (q) = vab {q)^bb (<?) - v 2 AB {q) [H]. Due 
to their generic form, the partial structure factors will diverge for all those conditions 
where D(q) = 0, and therefore the homogeneous fluid will be unstable. If the divergence 
takes place at q — 0, demixing occurs and the spinodal line denotes the locus of state 
points where -D(O) = 0. On the other hand, a divergence at some q* > signals an 
instability in the fluid (which is referred as Kirkwood instability) with respect to a 
periodic modulation of the density. The locus of state points where this instability 
takes place corresponds to the so-called A-line [15] . 

For the system at hand, the last scenario is indeed the case, where the instability 
dominating the phase behavior is mainly due to the Q ± -nature of the species B. Figured] 
displays in the (x, p)-plane the corresponding A-line estimated from MFA. Since this line 
is related to the limit of stability for the homogeneous fluid phase, it is expected that 
close to it the system will tend to crystallize into solids whose lattice constant will be 
dictated by q* (see Tabled]). By visual inspection of the configuration snapshots, state 
points for which the initial crystal structure was stable were discriminated from those 
for which the latter melted. The melting line estimated in this way was consistent with 
the analysis of Bragg peaks in the partial structure factors (see below). No transition 
to other crystalline structures was observed during the simulations, confirming the fee 
structure as the underlying lattice in the investigated crystalline states. The typical 
snapshots of Figure [1] indeed exhibit the expected ABC staking of the (lll)-planes. It 
should be noted that the apparent contradiction of having, in Figured] a fluid phase that 
is stable beyond the A-line is an artefact arising from applying the MFA in calculating 
the latter. The MFA becomes increasingly accurate at T grows, whereas we are working 
here at T = 0.3. Nevertheless, this has no consequences in what follows, since we are 
employing the MFA A-line merely as an indicator of the region in which alloy formation 
(crystallization) is expected, without basing any further quantitative predictions on its 
precise location. 
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Figure 2. Radial distribution functions (top) and structure factors (bottom) for 
x = 0.80 and different densities. The arrows in the plot of <?_BB(r) indicate the position 
of the fcth-neighbor in a fee lattice, i.e. dk = Vkd nn = \fkj2 af cc . The values of the 
selected densities are the same in all panels (see legend for §aa (r)). 



More detailed structural information can be gained by considering the partial radial 
distribution functions g a p{r\ which measure the probability to find a pair of particles, 
of species a and j3, at a mutual distance r. In figure |2] these quantities are displayed 
for x = 0.80 and different densities. As expected, QBB(j r ) exhibits a large peak at 
r = (namely gsBir ~^ 0) ~ 30), indicating the formation of clusters, whose size 
extends up to d c ~ 0.75. The second peak at r ~ 1.34 corresponds to the nearest- 
neighbor distance in the fee crystal, d nn = af cc / a/2. The radial distribution function for 
distinct pairs, gAB{r), displays rather different features. The correlation hole at r — > 
indicates that the small A-particles avoid the lattice sites and tend to localize in the 
interstitials between 5-clusters. Concerning correlations between pairs of A-particles, 
9AA{f) displays finite small values (~ 0.2) at r — > 0. This means that for the considered 
conditions, no appreciable clustering of the species A is induced by the clustering of the 
species B. 

On increasing the density, a significant structural change occurs in the system, 
which is evidenced by the emergence of well-defined peaks of all g a i3(r) at large distances 
(see top panels of figure [2]). This change is better reflected in the corresponding partial 
structure factors S a/ 3(q). As shown in bottom panels for the largest densities, the 
emergence of Bragg peaks in Sbb(q) indicates a rather ordered structure, which still 
retains some characteristics of the fluid phase. On the other hand, the nearly flat 
structure of Saa(q) clearly suggests the presence of a, nearly ideal, fluid phase of 
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A-particles. This fluid phase is immersed in the crystalline matrix of clusters of B- 
particles. Strictly speaking, the A-particles form a very weakly modulated fluid, so that 
the mixture of localized -B-particles and delocalized A-particles forms a sublattice melt 
phase. The one-particle density of the A-particles, however, cannot be strictly uniform, 
since the crystallized .B-particles act on the A-species as a periodically modulated 
external potential. This structural scenario has its dynamic counterpart, as will be 
shown below. 

The occupation number n of a given lattice site, i.e., the number of .B-particles 
in the associated cluster, is obviously an integer number. The average occupation is 
in general non-intee; 61* , £lS db result of (at least two) significant contributions of distinct 
values of n all over the lattice. For example, in the crystalline system of composition 
x = 0.80 and density p = 2.40 around the 65% and 25% of the clusters are composed 
respectively by three and four .B-particles. This percentages were estimated by counting 
the coordination number of each particle for r < d c ~ 0.75. On passing, it was also 
estimated by using the same criterion that around 90% of the A-particles are isolated 
from particles of the same species. The former results reflect that the lattice is full of 
defects which might be expected to break its stability. As discussed in Ref. p2] for 
the one-component cluster crystal, the stability of the lattice is actually maintained by 
incessant hopping of the -B-particles between distinct clusters. This is also the case in 
the mixtures considered in the present work, as discussed in the next subsection. 

3.2. Dynamics in real space 

A first step to gain some insight into the transport properties of the system is provided 
by the mean squared displacement (Ar^(t)) (MSD). This quantity is shown in figure [3] 
for different densities and relative compositions x = 0.65 and x = 0.95. In the case of 
the A-particles a rapid crossover is observed, for all the investigated densities, between 
short-time ballistic motion ((Ar^(t)) oc t 2 ) and long-time diffusion ((Ar^(i)) oc t). This 
is the usual dynamic behavior of a highly mobile fluid phase. Thus, as anticipated by 
the static correlations (see above), the mixture reaches a phase in which the species A 
shows a fluid-like behavior in a crystalline matrix of clusters of .B-particles. 

The MSD of the .B-particles shows a similar simple behavior in the fluid phase. 
However, an abrupt dynamic change is revealed by crossing the crystallization point. 
Thus, an intermediate plateau regime arises between the ballistic and diffusive limits, 
extending over longer time intervals as density increases. This regime is associated 
with the temporary trapping of the -B-particles within the clusters [12| EE] . The large 
oscillations at the beginning of the plateau regime are the signature of strong intracluster 
vibrational motion. Such short-time strong oscillations are an artifact of the simulated 
Newtonian dynamics. They are expected to be strongly damped under more realistic 
Brownian dynamics (BD); it has been shown that they even vanish in Monte Carlo 
(MC) simulations [IS] . Finally, at long times, the MSD displays a crossover to diffusive 
behavior. This scenario, which has been previously observed in the cluster crystal phase 




Figure 3. MSD of both species for x = 0.65 (left panels) and x = 0.95 (right panels). 
The arrows in bottom panels separate fluid states for the B-particles from crystal 
clusters (see also Figured]). 



of the pure GEM-8 system [T21 EE] , does not originate from the contribution to the MSD 
of rare events involving a few diffusing 5-particles. As anticipated above, it is the result 
of incessant hopping of 5-particles between neighboring lattice sites. Thus, the hopping 
mechanism fully changes the initial identity of each cluster without altering the overall 
crystalline lattice structure [12], [13] . 

This hopping mechanism is clearly reflected in the van Hove self-correlation, 
G$(r, t), of the .B-particles. Figure H] shows a representative example for composition 
x = 0.80 and density p = 3.05. A sequence of well-defined peaks is present in Gg\r, t), 
corresponding to different distances between lattice sites. With increasing time, the first 
peak in Gr (r, t) decreases while peaks located at larger distances grow progressively, 
corresponding to particles moving away from their original home clusters. By simple 
integration of Gg(r, t) from its first minimum to r — > oo, it is found that only the 30% 
of the -B-particles are located at t ~ 50000 in their initial home cluster. As expected 
from the observations in the MSD (see above) the van Hove self-function for the non- 
clustering species A displays a more simple behavior. For t ~ 50 the initial peak at 
r pd 0.3 has already vanished. At that time scale almost all the A-particles have left 
their initial interstitial positions and G^ (r, t) exhibits simple Gaussian behavior. 

We can estimate the diffusivity of each species from the relation D a = 
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Figure 4. Van Hove self-correlation function at different times for x = 0.80 and 
p = 3.05 (left: A-particles, right: B-particlcs). Distances are rescaled by d nn (see 
text). 



10" 



10" 



10" 



10" 



10" 



10" 



4_io r 



1.5 2.0 2.5 3.0 3.5 4.0 
, ^_ 




c^ x - 


0.65 


M X = 


0.80 


X = 


0.95 




1.00 














1 




8.0 



4.0 



5.0 



6.0 
px/r 



7.0 



9.0 



Figure 5. Diffusivities of the A- and -B-particles as a function of the density (inset) 
and of the combined variable px/T (main panel). Solid lines are guides for the 
eyes. The dashed line in the main panel indicates approximate Arrhenius behavior 
D B ~ cxp[-2.3/9i3/T] (see text). 



lim^oo (Ar^(t)) /6t. Figure [5] shows the obtained results as a function of the quantity 
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Pb/T = px/T. For each composition x, an abrupt drop in the diffusivity of the B- 
particles is observed in a narrow range of density. This drop reflects the transition of 
the species B from the fluid to the cluster crystal, and it is indeed expected from the 
previous observations in the MSD (see Figure [3]). The approximate scaling behavior 
for Db in the cluster crystal phase generalizes the observation for the pure GEM-8 
system [121 H3]- The latter was rationalized as the result of activated hopping motion 
between minima of the local potential energy, which are placed at the lattice sites 
[T2l IT3| [T9| 120] . It was found that the separating energy barrier between neighboring 
sites scales as AU ~ 2.3ps, leading to Arrhenius behavior D B ~ exp[— AU/T] [121 113] - 
Data of Figure |5] show that this observation is not altered in the mixture over a broad 
range of values of px/T. The introduction of the non-clustering A-particles does not 
even change significantly the former activation energy AU. 



3.3. Scattering functions 

Relaxation of density fluctuations of wave vector q are evaluated by means of the 
intermediate coherent and incoherent scattering functions. The coherent function is 
defined as 

p („a- (p a M)M-q,o)) (A . 

fap{q,t - 7 — / — pre — 7 tttt (4) 

(p a (q, 0)pp(-q, 0)) 

and characterizes collective a-(3 correlations, whereas the incoherent function accounts 
for self-correlations and is given by 

i /N a \ 
Fi s) (q,t) = — ^£exp{zq ■ [r QJ (t) - r aj (0)]}J • (5) 



With the used normalizations F a/3 (q,0) = F^(q,0) = 1. Figure [6] shows results 
of the scattering functions for composition x = 0.95 at densities above and below the 
crystallization line. The selected value of q = 4.0 is not a reciprocal lattice vector 
(RLV) and therefore data in this figure reflect relaxation of out-of-lattice correlations. 
Consistently with the observations in the real space (see above), both incoherent and 
coherent functions for the species A exhibit a fast decay at all the investigated densities, 
also for those where the S-particles form the cluster crystal. This is reminiscent of the 
scenario presented for the small particles in Yukawa mixtures immediately below the 
crystallization temperature [2T]. In close similarity to the aforementioned work, our 
data reveal, for the crystalline phase, a weak oscillation at t ~ 2, immediately after the 
microscopic decay. It reflects short-time vibrational dynamics of the A-particles within 
the interstitials. This effect was not detected in the respective MSD (see above), possibly 
because contrary to the scattering functions, the former is generally dominated by fast 
contributions, which may mask the oscillations. As discussed above, these short-time 
oscillations and the similar ones observed for the scattering functions of the 5-particles 
(Figures El and [7]) are artifacts of the Newtonian dynamics and will be strongly damped 
or even vanish for BD and MC simulations. 




Figure 6. Intermediate incoherent and coherent scattering functions for fixed x = 0.95 
and q = 4.0, at different densities (see legend for Fj± (q,t)). The arrows in bottom 
panels separate fluid states for the i?-particles from crystal clusters (see also FigureHJ. 



Concerning the dynamics of the clustering I?-particles, the present system exhibits 
more complex features than usual crystals. As in real-space quantities (see above), 
coherent and incoherent functions reflect an abrupt slowing down of the dynamics 
by crossing the crystallization point. The intermediate plateau regime extends over 
longer time intervals as density increases, before the ultimate relaxation at long times. 
This observation is reminiscent of the usual dynamic scenario associated to the glass 
transition. Having noted this, distinct features are revealed by inspection of incoherent 
and coherent functions. Thus, the relaxation of the coherent function Fbb{<1^) in the 
cluster crystal phase is consistent within statistics with a same plateau height (Debye- 
Waller factor, f q ) for all the investigated densities. This is not the case for the incoherent 
function F B s \q, t), where the plateau height (Lamb-Mossbauer factor, /W) clearly grows 
up as density increases. 

The usual observation in glass-forming systems is that both Debye- Waller and 
Lamb-Mossbauer factors remain constant on approaching the glass transition from the 
ergodic phase. On crossing the transition and moving deep into the glassy phase, the 
former factors f q and progressively increase. This feature reflects the progressive 
decrease of the localization length respectively associated to collective and self-motions. 
With these ideas in mind, the increase of the plateau height in the incoherent functions 
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for the -B-particles (Figure E]) suggests that a localization transition has occurred for the 
self-motions, tentatively at the crystallization density. On the contrary, such a transition 
is not detected for out-of-lattice collective correlations at any of the investigated 
crystalline states. This dynamic decoupling, in the meaning of a different locus for 
the localization transitions of self- and out-of-lattice collective motions, is also present 
in the pure cluster crystal [121 US]. Having noted this, the ergodicity of self-motions 
is, as in the pure GEM-8 system, restored at long times. This is reflected in the 
full decay of F#\q,t). In other words, the localization transition for self-motions of 
the -B-particles is actually avoided. The mechanism which restores ergodicity of self- 
motions is naturally the incessant hopping between neighboring clusters discussed above. 
Qualitatively similar findings are observed for the other studied compositions x = 0.8 
and x = 0.65 (not shown). Thus, the addition of non-clustering A-particles, even up 
to a 35% of the total, does not alter the former scenario of dynamic decoupling in the 
B-particles, and the (avoided by hopping) localization transition for self-motions still 
occurs at much lower densities that for out-of-lattice collective correlations!^ 

As discussed in Refs. [T2l [T3] , the former scenario has similarities and differences 
with the dynamics in plastic crystals, where molecules are constrained to vibrate around 
lattice sites, but perform full rotations which relax out-of-lattice collective correlations. 
In the case of cluster crystals, cluster deformation constitutes an additional mechanism 
of relaxation for such correlations. However, contrary to the case of plastic crystals, self- 
motions can explore arbitrarily long distances in cluster crystals by means of activated 
hopping. This leads to non-zero diffusivities (Figure [5]) and full relaxation of incoherent 
scattering functions (Figure [6]). 

As discussed above, the data presented in Figure [6] for a selected non-RLV q = 4.0 
display similar features, in the case of the species B, to those of the pure GEM-8 
system. Novel features are revealed for both species when the scattering functions are 
represented as a function of q. This is demonstrated in Figure [7] for a fixed crystalline 
state point, x = 0.95, p = 2.70, and for different values of q, including RLVs and non- 
RLVs. In the case of the A-particles, the incoherent function displays no signatures of 
glassy dynamics, but a fast decay for all the wave vectors. This is also the case for 
the coherent function at non-RLVs. Indeed the respective time scales are very similar 
for common q. However, a strong decoupling between self- and collective dynamics is 
observed for RLVs. Contrary to the incoherent case, the coherent function exhibits a 
plateau over several decades before relaxation at long times. 

This dynamic decoupling is rather different from that above discussed for the species 
B at non-RLVs. On the contrary, it is somewhat reminiscent of the scenario displayed 
by the small particles in disordered binary mixtures for certain ranges of composition 
and size ratio [22l [23] . In such systems the small particles move along a channel-like 
structure, formed by the interstitials of a matrix of large particles which relaxes in a 
much longer time scale, or which is even in the glassy state. As a consequence of the 

§ It should be mentioned, however, that although the ^4-particlcs constitute 35% of the number density 
of the system, they only account for less that 1% of its volume fraction due to their much smaller size. 
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Figure 7. Intermediate incoherent and coherent scattering functions for fixed x = 0.95 
and p = 2.70, at different wave vectors (see legend for F^'(q,t)). Solid lines and 
symbols correspond respectively to RLVs and to non-RLVs. The dashed lines indicate 
the zero value of the correlators. 

quasistatic arrangement of the interstitials, collective correlations of the small particles 
decay very slowly, despite fast self-motions are performed, allowing for the exploration 
of large distances. Thus, decoupling of incoherent and coherent scattering functions of 
the small particles is observed in all the range of low and moderate wave vectors which 
probe the matrix structure |22j [23] . As discussed above, for the species A of the mixtures 
here investigated, the former decoupling is only observed for RLVs. This reflects the fact 
that the non-clustering particles are preferentially located at the interstitials between 
the .B-clusters, as was seen in the static correlation functions (Figure [2]), and therefore 
also move preferentially between them. Since collective correlations between interstitial 
positions are probed by the RLVs, a decoupling between incoherent and coherent 
dynamics is observed for such wave vectors, in analogy with the scenario discussed 
above for disordered mixtures. 

Concerning the scattering functions of the clustering species B, slow relaxation is 
observed for the incoherent and coherent cases, and both for RLVs and non-RLVs. As 
usual, Lamb-Mossbauer factors show a monotonic g-dependence, while Debye- Waller 
factors follow the modulations of the structure factor. Having noted this, coherent 
functions reveal non-trivial differences between RLVs and non-RLVs at the longest times 
of the simulation window. Namely, FBB(q,t) exhibits, for RLVs, an oscillation around 
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Figure 8. Intermediate coherent scattering functions of B-particles for several RLVs 
(see legend). Filled symbols correspond to the mixture with x = 0.95, p = 2.70, and 
T = 0.30. Empty symbols correspond to the pure GEM-8 system with p = 3.0 and 
T = 0.40. Solid lines are guides for the eyes. The dashed line indicates the level 
F BB (q,t) =0. 

t ~ 10 5 . This oscillation is apparently absent for non-RLVs. 

For the sake of clarity, the former data for the RLVs are shown separately in Figure [8] 
over a time window t > 10. Given the large amplitude of the long-time anticorrelations, 
reaching values of Fbb{q, t) = —0.2 at the minimum, it is improbable that the oscillation 
around t ~ 10 5 is an statistical artifact. The figure also includes coherent data, at its 
respective RLVs, of the pure GEM-8 system at p = 3.0 and T = 0.40. The rather 
different behavior in the pure system and in the mixture, even with only a 5% of A- 
particles in the latter, becomes evident. While full relaxation followed by oscillation 
is observed in the mixture, no decay is present in the pure cluster crystal. It must be 
stressed that this difference is not related to a much slower intrinsic dynamics of the 
pure system for the selected state point p = 3.0 and T = 0.40. Indeed, by inspection 
of all the other dynamic observables introduced above, including coherent functions for 
non-RLVs, it is found that the pure system shows a faster dynamics than the mixture 
at the former state points. For example Db ~ 2 x 10~ 6 and 10 -5 respectively for the 
mixture and the pure system. Therefore, the oscillatory feature observed at long times 
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in the mixture is induced by the addition of the gaussian A-particles. 

The presence of negative correlations in the collective dynamics of S-particles at 
near RLV-wavevectors, reached at typical times r a ~ 10 5 , is a particular feature of the 
mixture, which is is absent for the pure B system within the simulated time window. 
Physically, it describes a process by which a positive density correlation at t — at 
a distance of the order of the lattice constant turns negative at time scales r a . As 
such, it points to the existence of some characteristic oscillation frequency u a ~ r~ l , 
which corresponds to wavevectors lying at the edge of the Brillouin zone, i.e., a short- 
wavelength acoustical mode. Although such a process has not been seen for the pure 
-B-system, phonon excitations exist for the pure system as well; it is thus reasonable to 
assume that the corresponding frequency for the pure system is much smaller that u a , 
so that the time required for this relaxation process to be seen in the coherent scattering 
intensities is much longer. A detailed calculation of the phonon spectra of the mixture 
at hand, along parallel lines to the one recently carried out for pure systems [21] would 
shed light into this question, it is however cumbersome and it will be left as the subject 
of a future publication. Nevertheless, in view of the fact that the A-particles provide 
additional repulsions to the £>-species and thus enhance the restoring forces that act on 
the latter, it appears plausible that in the mixture the phonon frequencies are larger 
than those in the pure system and thus the characteristic time r a becomes visible within 
the simulation window. Alternatively, the issue could also be resolved by performing 
very long simulations for the pure S-system. 

4. Conclusions 

The dynamic aspects of binary mixtures of ultrasoft (GEM-8 and gaussian) particles 
have been studied by means of MD. The present work extends previous results for 
the cluster crystal phase of the pure GEM-8 system, and investigates the effects of the 
addition of the non-clustering gaussian particles on the corresponding dynamic scenario. 
The obtained results show that as the total density increases at fixed composition, the 
GEM-8 species builds cluster crystals with a relatively high localization. As in the one- 
component system, the incoherent scattering functions indicate a localization transition 
for self-motions of the GEM-8 particles, which is avoided by incessant hopping motion 
between clusters. This transition seems to occur at lower densities than for out-of-lattice 
collective correlations, confirming in the mixture the dynamic decoupling observed in 
the one-component cluster crystal. 

On the other hand, the gaussian particles remain delocalized in the periodic 
potential induced by the GEM-8 clusters and display fast self-motions. However, slow 
collective dynamics is observed for specific wave vectors, namely those belonging to the 
reciprocal lattice. This feature reflects the preferential motion of the gaussian particles 
over the interstitials, confirming the expectations from static correlations. A striking 
feature is revealed by the analysis of collective correlations of the GEM-8 particles, for 
wave vectors at the reciprocal lattice, a feature which we attribute to the stiffening of 
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an acoustical mode at the edge of the Brillouin zone, caused by the presence of the A- 
particles. An additional open question that should be the subject of future investigations 
is the influence of the type of elementary processes employed to model the dynamics in 
the mixture, along the lines of the work by Coslovich et al. carried out for pure systems 
[18]. Here, it should be examined which of the features discovered by means of our 
MD approach survive if one employs Monte Carlo moves instead, which are expected to 
mimic better the overdamped, Brownian Dynamics of the system. 

Finally, while there is evidence that both classes of ultrasoft potentials accurately 
model the interactions between real macromolecules in dilute systems, the validity of 
these potentials in dense systems has yet to be verified. The results presented in this 
paper add to a growing body of work which suggests that, the search for dense systems 
where ultrasoft potentials remain valid is a worthwhile pursuit. 
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